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Abstract 

The D2-D0 bound state exhibits a Gregory-Laflamme instability when it is sufficiently 
non-extremal. If there are no DO-branes, the requisite non-extremality is finite. When most 
of the extremal mass comes from DO-branes, the requisite non-extremality is very small. 
The location of the threshold for the instability is determined using a local thermodynamic 
analysis which is then checked against a numerical analysis of the linearized equations of 
motion. The thermodynamic analysis reveals an instability of non-commutative field theory 
at finite temperature, which may occur only at very long wavelengths as the decoupling limit 
is approached. 
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1 Introduction 

The tension of the D2-D0 bound state is 



TD2-D0 — 




Here tdi and tdo are the D2-brane tension and DO-brane mass, respectively. N 2 is the number 
of D2-branes, and po = N0/V2 is the number of DO-branes per unit world- volume area. We are 
considering the bound state at extremality in asymptotically flat ten-dimensional spacetime. 

The D2-D0 bound state is supersymmetric and therefore stable. Without referring to 
super symmetry, we can see from (JTJ) that the bound state is stable against infrared pertur- 
bations that make the number density of DO-branes non-uniform: the key fact is that TD2-D0 
is a convex function of p . More precisely, suppose we make a long- wavelength perturbation, 
Po —> Po + Sp . Then the change in the total mass of the bound state is 

6M^ J d 2 i y/g + (Po + 5po) 2 r£ - ^V?r£ 2 + p^o) • (2) 

The right hand side is always positive, provided / d 2 a^fg5po = 0, precisely because the 
tension formula (JTJ is convex in p 2 . The equality in (J2J is approximate because there should 
also be contributions expressed in terms of spatial derivatives of 5po. The point of considering 
only infrared fluctuations is to suppress these contributions. Within this framework, we 
have reproduced the expected conclusion that the D2-D0 bound state is stable against local 
rearrangements of DO-brane charge. 

The purpose of this paper is essentially to extend the above analysis to a non-extremal 
D2-D0 bound state, where it is not so obvious whether there should be an instability. All 



1 



calculations will be done using the supergravity description of the bound state, so the quan- 
titative results are limited to supergravity' s regime of validity However, I expect that the 
following qualitative picture holds more generally for D2-D0 bound states without angular 
momentum: 

• Far enough from extremality, there is an instability. Except in special limits, the 
instability disappears when the mass above extremality is some 0(1) multiple of the 
extremal mass. 

• When DO-branes make up most of the mass of the bound state, the instability persists 
close to extremality: the aforementioned 0(1) factor becomes small. 

One of our main tools is a notion of local thermodynamic stability. Briefly, at finite temper- 
ature, we may independently rearrange DO-brane charge and entropy, so to redo the above 
calculation properly, we should demand that the energy is a convex function of them both 
in order to avoid an instability. More precisely, the Hessian matrix of second derivatives of 
the energy with respect to the DO-brane charge and the entropy should be positive definite 
to avoid an instability. 

In m |2] it was conjectured that for systems with translational symmetry and infinite 
extent, a Gregory-Laflamme (GL) instability [SJ E] arises precisely when the system has a 
thermodynamic instability — meaning that the Hessian matrix of second derivatives of the 
mass with respect to entropy and any charges which are capable of being redistributed over 
the directions in which the GL instability is supposed to occur. I will refer to this as the 
correlated stability conjecture (CSC). A heuristic justification of the CSC is that if one writes 
entropy in terms of the mass and charges, then a positive eigenvalue of the Hessian matrix 
identifies a direction of evolution in which entropy can increase while total mass and charge 
remain constant. The idea that the GL instability is associated with entropy increase was 
already present in the original papers pQUj. The CSC makes this idea more systematic by 
equating a local thermodynamic instability to a perturbative dynamical instability. 

For classical systems with horizons of finite extent (e.g. a Schwarzschild black hole, or 
a black string wrapped on a circle), there can be no conclusive argument for the existence 
of a horizon instability based solely on thermodynamics. The reason is not that infinite 
volume is necessary for a thermodynamic limit to be taken: with the freedom to rescale the 
Planck scale, any classical horizon can be regarded as large. Instead, the reason is that the 
wavelength for the instability may be larger than the system size. The generic expectation 
(with interesting exceptions near the boundary of stability) is that the wavelength of the GL 
instability is of the same order as the horizon radius. The claim, then, for finite size systems, 
is that the existence or non-existence of a GL instability is driven entirely by the competition 
between thermodynamic and finite-size effects. From this point of view, the Schwarzschild 
black hole in four-dimensions is stable only because of 0(1) effects: it's just a bit too small 
for its natural GL instability to fit on its horizon. 

A reason to be particularly interested in the GL instability for the near-extremal D2-D0 
bound state is that it interpolates between near-extremal D2-branes, for which it is known 
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there is no instability j3] and near-extremal smeared DO-branes, for which it has been argued 
that there is a GL instability @[Z1IB]. The interpolation itself is interesting: the dynamics 
of the bound state has been argued to correspond to non-commutative field theory (NCFT) 
jH]. More precisely, a particular limit, in which the energy scale of allowed processes is 
lowered at the same time as the density of DO-branes is increased, results in a decoupling 
of closed strings and excited open strings from the non-commutative dynamics of massless 
open strings. 

Because the NCFT limit corresponds to an arbitrarily large density of DO-branes, the 
convexity property that "holds together" the bound state is extremely weak: inspecting 
(JIJ for large po, we see that it is barely convex at large po- So we might think that a 
very slight non-extremality causes an instability. The question is whether the non-extremal 
temperature is bigger or smaller than the energy scale of allowed excitations in the NCFT 
limit. I will argue that it is smaller. Should we then conclude that NCFT is unstable at any 
non-zero temperature? The instability corresponds to the non-commutativity parameter 
becoming inhomogeneous. It's possible that the wavelength of this instability diverges in the 
NCFT limit. In other words, NCFT might avoid the instability by pushing it to infinite 
wavelength. This scenario seems likely based on decoupling arguments, but I do not have 
an explicit computation of the wavelength. 

In section El I briefly present the supergravity background describing the D2-D0 bound 
state. In section El I use the CSC to predict for which range of parameters a GL instability 
should occur. In section 0] I argue that the transition temperature to a GL instability 
vanishes as one approaches the NCFT limit. In section El I perform a numerical analysis 
on the equations of motion linearized around the D2-D0 background which confirms the 
predictions of the CSC. Section El comprises some concluding remarks. 



2 Classical action, solutions, and thermodynamics 



The part of the classical type IIA supergravity action relevant for the considerations of this 
paper is 



S = ^ f d l0 x VG 
Ik 1 J 



-20 



R + 4(d<f>)' 



-Hi 



If" - Ifi 



(3) 



with F A = F 4 + Ai A H 3 
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The non-extremal D2-D0 bound state corresponds to the following supergravity background: 



ds 2 str = H-^i-hdt 2 + D{dx\ + dx 2 )) + H 1 ' 2 [j-dr 2 + r 2 dfll 



A\ = coth a sin 0(1 — — ) dt 



1 

h l 

\ d± 

H, 

A 3 = coth a sec 9 (l - ^\ dt A dx 1 A dx 2 (4) 
e 20 = H i/2 D B2 = tan d f 1 __^ dx i A dx 2 

H ^ l + rtsmh 2 a D= 1 fc = i _ d . 

r 5 if 1 sin 6* + cos" 1 p r 5 

The thermodynamic quantities of interest are 



M=^r 5 (6 + 5sinh 2 a) 



47rr cosh a k 1 (5) 

/i2 = \i cos ^ Q2 = Q cos 6 1 /i = A 4 sin 9 Qo = Q sin 6* 

A = tmhQ 0= ^ sinhtt eo S ha, 

Ik 

where fi 6 = 7r 3 is the volume of a unit S* 6 . 

Solutions of the type described in (jlj are well known: see for example [TTH ITT] for sub- 
stantial generalizations of it, for references to earlier literature, and for some description of 
how such solutions can be obtained from simpler ones using rotations and T-duality. 



3 Consequences of the correlated stability conjecture 

To apply the correlated stability conjecture (CSC), we need to know what extensive thermo- 
dynamic quantities are capable of being locally redistributed. For the non-extremal D2-D0 
bound state, the D2-brane charge cannot be redistributed because D2-branes can't break. 
But the DO-brane charge can be redistributed. Also, entropy can be redistributed by making 
the horizon wavy. So to determine local thermodynamic stability, one should consider the 
Hessian of M with respect to S and Qo, holding Q2 fixed. 

In (jlj), M is not expressed as a function of S, Q2, and Qo] instead, these four quantities 
are expressed in terms of r , a, and 9. So it's useful to recall a fact from multi-variable 
calculus: if there is a smooth, invertible relationship between n variables Qi and n other 
variables qi, and M is known as a smooth function of the qi, then 

(dM_\ = d{Q l ,...,Q l ,M,...,Q n )/d{g l ,...,q n ) 
\dQi) Qj d(Q h Q n )/d(q u ...,q n ) 
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det H = O _ 9o2 i2t/2/c ^ ~ u o \ (~ ^ ~~ ^ cos 20 + 6 cos 2 6 1 cosh 2a) . 



where the denominator is the Jacobian, det(dQi/dqj), and the hat notation in the numerator 
is meant to indicate replacing Qi with M. Second derivatives of M can be taken by plugging 
a first derivative back into (JHJ) in place of M. Actually, the thermodynamic dual quantities 
T, /i 2 , and fio should be precisely the first derivatives of M with respect to S, Q 2 , and Qo 
because of the first law of thermodynamics: 

dE = TdS + fi 2 dQ 2 + fi dQ . (7) 

This consistency check was already made in ^Uj. The Hessian may be expressed as H = 
d(T, ji )/d(S, Qo), and the right hand side may be evaluated by repeated use of (JBJ). One 
finds 

k 4 sech 4 a 
87T 2 n 2 6 rl 2 V 2 2 (5 + 7cosh 2a) 

The boundary of local thermodynamic stability is where det H = 0. (More properly, it is the 
boundary of the region where all eigenvalues of H are positive, but it turns out in this case 
that H has at most one negative eigenvalue). Starting from JSJ) one can easily demonstrate 
that det H = is equivalent to 

csch a = \/3 cos 6 . (9) 

This marginal stability condition can be recast in terms of the potentials /io and fi 2 appearing 
in (jSJ) as 

^ + 4^=1. (10) 

To understand this result, it is useful to plot Qq/M = 5 sin 9 sinh a cosh a/ (6 + 5sinh 2 a) 
against Q 2 /M = 5 cos 9 sinh a cosh a/ (6 + 5 sinh 2 a), which is done in figure ^ Given the 
curve (jHJ) of local thermodynamic stability, the CSC predicts a GL instability precisely when 
a < csch -1 (-\/3 cos 9). If there are only DO-branes and no D2-branes (that is, 9 = 7r/2), then 
there is supposed to be a GL instability for all finite values of a. 



4 The non-commutative field theory limit 

Special interest attaches to the low-energy dynamics of massless open strings on the D2-D0 
bound state in the limit where most of the mass of the bound state comes from the DO-branes. 
This dynamics has been argued to be a non-commutative field theory (NCFT) decoupled 
from gravity and excited string states[S]. Yet we can see from figure Q that the boundary 
of stability predicted by the CSC goes through precisely the region of interest for this limit: 
Qo/M fa 1 while Q 2 /M <^ 1. The GL instability predicted by the CSC involves redistributing 
the DO-branes, which means making the parameter $ that measures non-commutativity a 
function of space — and of time unless there is an endpoint to the dynamical evolution of the 
GL instability. Such a situation would seem to spoil in an important way the decoupling 
of the NCFT from closed strings: spatially variable B 2 is clearly a closed string degree of 
freedom. 
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Figure 1: The thick curve indicates the boundary of stability. The shaded region is the 
unstable region. The thin curves show constant M for fixed Q = ^jQl-i- Qo'- the outermost 
of these curves is at M = Q, the next is at M = 1.15Q, and so on with even spacing 
continuing to M — 2.5Q. The vertical axis and horizontal axis show, respectively, the 
fraction of the mass that comes from DO-branes and D2-branes; the rest of the mass is 
non-extremality. 

How might the NCFT limit avoid the GL instability? One possibility hinges on the fact 
that it is a low-energy limit. The boundary of stability shown in figure ^ defines a critical 
temperature T c at any fixed values of Q 2 and Q . At least near the upper left corner of 
figure H when DO-branes make the dominant contribution to the mass, solutions on the 
stable side of the curve have temperature T < T c . If T c were large compared to the energies 
permitted in the NCFT limit, then clearly the GL instability would be altogether avoided in 
the NCFT limit. Unfortunately, this is not what happens: as we will show in this section, 
T c is small compared to the intrinsic energy scale of the NCFT, where i? is the non- 

commutativity parameter. 

There is another possibility: the critical wavelength above which there are unstable modes 
might diverge in the NCFT limit, for all temperatures which remain finite as compared to 
the scale 1 / \f&. In principle we should be able to check this numerically, but we have not 
done so. In section 15.31 we will briefly describe why this is a difficult numerical problem in 
practice. Without hard evidence one way or another, the divergent wavelength possibility 
is the one I regard as the most likely. If true, it implies that if we compactify the D2-D0 
bound state on a torus whose size is held fixed as one takes the NCFT limit, then there is 
no GL instability in this limit at all. A third possibility — that NCFT is unstable toward 
non-uniform non-commutativity parameter at some finite temperature — would simply mean 
that the NCFT limit is ill-defined. I would find this hard to understand in light of the very 
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plausible decoupling arguments put forward in 0. 

Let us now show that the critical temperature T c is much less than the intrinsic scale 
\j\p$ in the NCFT limit. The demonstration mainly consists of unraveling the definitions 
for T c and •& and recalling how the NCFT limit is taken. 

Following P, consider the D2-D0 bound state in the test brane approximation. The 
metric felt by closed strings is the metric Gmn appearing in (|2J). In discussions of NCFT 
this closed string metric is usually denoted Qmn^ while Gmn is reserved for the metric felt 
by open strings. I will follow this convention in this section. A useful quantity is 

E a p — 9a/3 + B a p , (11) 

where indices a, j3 denote a pull-back to the 2 + 1-dimensional brane world- volume. The 
inverse matrix E a/3 may be expanded into a symmetric and anti-symmetric part: 

E a P = G a P + — . (12) 
2na' 

Then G a p (defined as the inverse of G a/3 ) is indeed the metric felt by the open strings: this 
means, for instance, that massless open strings follow null geodesies of G a p. It is convenient 
to change coordinates on the brane to x a such that G a p = diag{ — 1, 1,1}. One finds that 
measures non-commutativity: 

[x a ,x p ] = i$ a(3 . (13) 

Operationally, this means that open string scattering amplitudes come with phases like 
e -^ al3 p a qi3 w here p and q are external momenta — see [§] for details. Assuming only space- 
space non-commutativity, the NCFT limit is a' — > with d = d 12 held fixed. Then \fd is 
an intrinsic length scale characterizing the NCFT. 

Consider the case where only _B 12 = — B 2 i = — tan 9 is non-zero. 1 One can easily show 
that t = t, x 1 = x 1 1 cos 9, and x 2 = x 2 / cos 9. Note that this means that the speed of 
propagation for massless open strings is cos 9. When expressed in the x coordinate system, 
$ 12 = — $ 21 = $ : where 

$ = 2Wtan#. (14) 

In the NCFT limit, i9 is held fixed while a' — * 0, so one needs 9 — > tt/2. 

The parameters used to specify the NCFT are the string coupling g s , the number iV 
of D2-branes, and the temperature T. Our aim now is to express the condition JSJ) defining 
the boundary of stability in terms of these quantities. First note that Q/V2 is the tension at 

1 The solution as expressed in nas -S12 = far from the brane and _B 12 f« tan0 at the horizon, where 
the approximate equality becomes exact in the near-extremal limit. A global shift, B12 — » B\i — £?, is 
appropriate to compare with the test brane approximation. 
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extremality, so Q 2 /V 2 is the tension at extremality of the D2-branes only: 

— = -z-^ r o smh a cosh cos u = N 2 T D2 = — , 9 /C!/9 • (15) 
V 2 2 k g s 4:T[ z a' 6 i z 

Using this equation together with the standard relation 2k 2 = (2ir) 7 g 2 a' 4 , one may eliminate 
ro in favor of g s , N 2 , a, 8, and a'. Using in addition (fT^j) . one may express the dimensionless 
quantity \fdT in terms of the same parameters: 

^ = 217/1 J /%9/10 { J 2)US ^^ (coserV^Te. (i6) 

Note that T is an energy scale measured with respect to t = t, while \fd is a length measured 
with respect to x 1 and x 2 . The correct dimensionless quantity to compute is \f&T /v, where v 
is the propagation speed of massless open strings. But we are using the coordinates (t, x 1 , x 2 ), 
both to define $ (see (O) and to measure T. And in these coordinates, v = 1. So \fQT 
really is the right dimensionless quantity after all. 

Now let us find the value of \f&T corresponding to the boundary of stability: using (jH)) 
to eliminate a from (J16)) one obtains 



JtfT - 5 ' 31 10 - fcot^ 3/1 ° 

V ^ C " 217/10^9/10 (^1/5^^ 



1 - -cot 2 fl + 0(cot 4 i 



^1 i/wy/r 

217/10^9/10 (^AT 2 )l/5 \ # I 



(17) 

l + 0(a a /^) . 



As described near the beginning of this section, T c is a transition temperature, below which 
there is a uniform phase, and above which there is a GL instability which makes $ non- 
uniform and whose endpoint is not known. 

From (fTTj) it is clear that \f&T c — » in the limit where a 1 — > 0, assuming finite 't Hooft 
coupling g s N 2 . This seems to invalidate the claim fT2\ IT3] that NCFT has the same large 
N thermodynamics as the corresponding commutative theory: introducing an instability for 
all T > is a pretty radical change! But if the divergent wavelengths hypothesis described 
near the beginning of this section is correct, the apparent discrepancy is merely an order of 
limits issue: the claims of ITB*] would follow if we take the NCFT limit on a torus of fixed 
size and then make the torus large, whereas to keep the instability in the theory one must 
keep the torus larger than the critical wavelength as one approaches the NCFT limit. 



5 Perturbing the D2-D0 bound state 

Although well-motivated and well-checked for a wide class of examples, the CSC is not proven 
in full generality. I therefore aim to demonstrate explicitly that when the CSC says there 
should be a GL instability for the D2-D0 bound state, there is a stationary perturbation 



S 



of the supergravity solution (j3J) that is non- uniform in the coordinate x 1 . It is generally 
understood (see for example ^3]) that if such a perturbation exists at a finite wavelength, 
then similar perturbations with longer wavelengths will be unstable. 

Obtaining the explicit form of the stationary perturbation requires numerics. The super- 
gravity action and background, © and (j3J), are somewhat complicated, and the perturbation 
equations are remarkably complicated. Therefore I will only present a summary discussion 
in this paper. Briefly, in section 15.11 the problem is reduced to two dimensions; then in sec- 
tion l5.2l it is explained how one extracts from the two-dimensional lagrangian and constraints 
a well-posed initial value problem. Finally, in section 15.31 1 describe the results of numerics. 

5.1 Reducing to a two-dimensional lagrangian 

I will assume that this perturbation is an s-wave with respect to the S 6 : it depends only on 
x l and r. Therefore it is natural to start by making a Kaluza-Klein reduction on the S G , 
t, and x 2 directions. The result is an effective lagrangian for several scalars, equation (|29|) . 
together with constraints from the two-dimensional Einstein equations described in equation 

The ansatz may be described as follows: 

ds 2 str = G ilv dx^dx v - e 2Lpl dt 2 + e 2{p2 dx 2 2 + e 2 ^dVt\ 
B 2 = b 2 + bi A dx 2 + b 1 Adt + b dt A dx 2 

H 3 = h 2 A dx 2 + h 2 A dt + hi A dt A dx 2 
Ai — a\ + aodx 2 + dodt (18) 

F 2 = f 2 + fiA dx 2 + h A dt 
A 3 = c 2 A dx 2 + c 2 A dt + ci A dt A dx 2 

F 4 = f 2 Adt A dx 2 , 

where h 2 = dbi, h 2 = db±, and so forth. Forms such as h 2 and bi are constructed from the 
dx a , where x a runs over (x l ,r). 

Plugging (J18)) into the ten-dimensional action in (@J), one obtains a two-dimensional action 
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that is the sum of two terms: 



(19) 



S = J d 2 x(L NS + L} 
G- l ' 2 L NS = e~ 2<1> (r + 4(«9$) 2 - (d Vl ) 2 - (d<p 2 ) 2 - 6(c^ 3 ) 2 + 30e" 2 ^ 

_ l - e - 2 ^h 2 2 + l -e~ 2 ^hl + ie" 2 ^" 2 ^ 

it ^ 2 2 

+ i e -^-^+ 6 ^(/ 2 _ ao /, 2 + a /i 2 + ai A . 

Here G = — det G Q/ 3 and R refer to the two-dimensional string metric, and the two-dimensional 
dilaton is 

$ = 0--(^ 1 + ^ 2 + 6^ 3 )- (20) 

The Wess-Zumino term in the action in (jlj) vanishes for the ansatz (J18|) . 

The fields in (fT9*|) include four gauge fields: b±, b\, a±, and ci. These fields are non- 
dynamical: they can be eliminated in favor of constants of the motion which are roughly 
their conjugate variables. To illustrate this process of elimination, consider the following 
simple example: 

S = J C , C = ]f{h + Q2) A *(/ 2 + q 2 ) + g / 2 , (21) 

where ip, g 2 , and go depend arbitrarily on other fields. The manipulations for eliminating f 2 
are as follows: 

5£ = e^dSai A *(/ 2 + q 2 ) + q^dSai = ir a d5ai 

7r a = e v * (f 2 + g 2 ) + g (22) 
C = C- 7i a dai = -^e" v * (7r a - g ) 2 + (n a - g )g 2 . 

The "reduced lagrangian density" £ is the analog of the Routhian in classical mechanics. If 
one adds to the original action terms independent of a\ and its derivatives, but depending 
on other fields which go and g 2 may also depend on, then the equations of motion for these 
other fields may be obtained by varying the reduced lagrangian. 
Starting with the lagrangian density for the form fields, 

C f = --e^^h, A *h 2 + - e ~ 2 ^ 2 ^h 2 A *h 2 + l e -™-*n-*P*h x A *h x 
2 2 2 

_ _ e ^i+^2+6^ 3 j 2 A ^ _ 1^1-^2+6^3^ A ^ + Ig-Vi-Hw-Wvs^ A # £ (23) 

+ I e -^-^+ 6 ^(j 2 _ ^ + a Q h 2 + ai A hi) A *(/ 2 - a /i2 + d h 2 + a a A hi) , 
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One can eliminate first Ci, then a 1; then b\ and bi in favor of the following constants: 

-2*-2C£?2 !, ~ ~ ^ ' 

7r fe = - * e ^ 2 /i 2 + vr c a 
7r b = *e ^ ii 2 — ir c a . 

The reduced lagrangian density is 
t = * (tt, - 7r c a ) 2 - l -e 2 " +2 ^ * (tt 6 + ^ c a ) 2 + I e -<«-<»-<*» * - ^ ) 2 

1 _L,„ fi,„ O 1 _0(T>_0,n. r 7 1 



-6<^3^.2 



0; 



e + -e- 2 *"^-^^! A *~h - - e ri-w+s<ps fl A +/l (25) 
+ ^e-^+^A A *A . 
The total reduced lagrangian (restoring the first line of L^s in (JISJ) is 

G-V2 L = e -2* ^ + 4((9$) 2 _ ( ^ i)2 _ ( ^2 _ 6( ^ s) 2 + 1 f 

- ^-^^(dao) 2 + ^e^ 1+ ^ +6 ^(<% ) 2 

+ We- 2 *-** + V +2 ^K - 7f c a ) 2 - V+ 2 ^ + ^r c a ) 2 
+ -e-^i-^-e^^ + ^ Q )2 _ 1 ^+^-6^3^.2 _ 

Clearly, this lagrangian describes two-dimensional dilaton gravity coupled to the seven scalars 
($, ipi, ip 2 , <P3, bo, a , a ), all of which participate in the potential energy function. 

From (J3J) and (|18|) one can easily read off the fields involved in the background solution: 



(26) 



dsn tf- = —^dx-i H — dr 

' JH h 



e 2 ^ = ^= e 2 ^ = -S= e 2 ^ = VHr 2 e 2 * 



H VH \hHr & (27) 

ao = cothasin^ ^1 — — ^ ci = cothasec6> ^— — 1^ dx 

6i = tan6» ^1 - — ^ dx 1 . 



The other fields in the action (JT9)l vanish. One easily obtains 

7f c = — 5?"q cosh a sinh a cos 6* 7r a = 
7Tb = — 5tq sinh 2 a cos 6 sin # 7Tfc = . 



(2* 
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It is straightforward though somewhat tedious to show that (|27p satisfies the equations of 
motion derived from (|26|) . 

In obtaining the linearized perturbation equations, it is efficient to introduce a conformal 
factor: send G a p — > e 2a G a p, so that the total reduced lagrangian becomes 



(db ) 



G-^L = e~ 2 * [r - Adad® + 4(<9$) 2 - (9^) 2 - (d<p 2 ) 2 - Q(d<p 3 ) 2 + X -e 

- -e 2a+Lpl ^ 2+ ^(da ) 2 + -e 2a -* 1+ * 2+ ^{da ) 2 
2 2 

+ ^-^-^ + i e 2CT+2 * +2 ^(^ - ^o) 2 - ^e M+2 ^(^ + n c a ) 2 

+ I e 2<T - Vl -^- 6v3 (vr a + ^ c 6 ) 2 - l e 2 -+^i+^-6^3^2 _ 
2 2 

(29) 

Eight second order equations of motion are obtained by varying with respect to the eight 
scalars a, $, <pi, <p 2 , f3, b , a , and a . It is not helpful to write out these equations in detail 
here. One also has the Einstein equations, obtained from varying with respect to the metric. 
If we compress f!29|) to the notation 

G -i/2 L = e -2* (jR + 4(9$)2) _ \g abdrd(j)b _ ^ (30) 
where 4> a includes all eight scalars, then the Einstein equations read 



'~ 2 * ( iV - \RG^ + 2G^(d^) 2 + 2V At ^$ - 2G^V 2 $ 



(31) 



The Einstein tensor R^ — ^RG^ U vanishes because we're working in two dimensions. One 
combination of the Einstein equations (the trace) is equivalent to the equation of motion 
one gets by varying a. The other two independent combinations can be viewed as gauge 
conditions for the gauge choice of writing the perturbed metric as a conformal factor times 
the original one. 



5.2 Setting up a well-defined initial value problem 

Because our interest is in a stationary mode with some wave-number k, let's restrict the 
ansatz by letting each of the eight scalars take the form 

(ft = 0o ( r ) + S(j)(r) cos kx , (32) 

where is one of the eight scalars and 0o( r ) is the background solution. Terms proportional 
to sin kx turn out to be unnecessary: they decouple from the cos kx perturbations at first 
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order. Plugging (|32|) into the equations of motion and constraints and expanding to linear 
order in the S(p(r) leads to linear ordinary differential equations. These linearized equa- 
tions comprise eight order equations of motion (obtained by varying the two-dimensional 
lagrangian (|29j)) and two constraints (obtained from the Einstein equation) in the eight vari- 
ables a, $, tp%, <f2, <fi3, bo, do, and a . The form of these ten equations is still remarkably 
complicated and will not be reproduced here. 

The ten equations described incorporate some redundancy: six of the equations of motion 
plus the two constraints imply the remaining two equations of motion. One can therefore 
drop the equations of motion for 5a and 5$. It is found by direct computation that the 
equations of motion for Sao and 5bo decouple from the others, and that these two fields may 
be set uniformly to zero. What remains is six coupled equations in the six variables 

qi = (5a, <5$, Sift, 6<f 2 , S<fi3, 5a ) . (33) 

They comprise the equations of motion for 5(fi, 5if2, 5ip3, and 5ao, plus the two constraints 
from the Einstein equations. These linear equations have terms multiplying qi(r), q[(r), 
and q"(r), except that 5a"(r) and 5<3>"(r) do not appear. Let us call these six differential 
equations the perturbation equations. 

Boundary conditions at infinity are simple: one should require that the qi are normaliz- 
able. For k ^ 0, the behavior for large r is q^ ~ e likr for some constants 7$. A normalizable 
solution is one in which all the 7^ are negative. 

Boundary conditions at the horizon are also conceptually simple: one should require 
that the perturbed metric and matter fields are smooth at the horizon in ten dimensions. 
This requirement becomes rather obscure in our two-dimensional language. So I implement 
boundary conditions in a different way, which is somewhat complicated to state, but which 
I expect is equivalent. Start with the eight non-trivial equations obtained by setting 5a = 
5b = in the original ten equations described below (}3*2*|) . Then set r = tq. It happens 
that all terms involving qf(ro) drop out: the second derivatives are multiplied by functions 
which vanish at r = r . So the result is eight linear algebraic relations on the twelve 
quantities (?i(^o); ^( r o))- Boundary conditions at r = r on differential equations which do 
not involve 8a" (r) and 5$"(r) should not themselves depend on §a'(ro) or 5$'(ro). Imposing 
this limitation, one goes from eight relations to five, now involving the ten quantities 

5a, 5$, 5(fi, 5(p 2 , 5ip 3 , 5a , 5ip[, 5ip' 2 , 5ip' 3 , 5a' . (34) 

Of these conditions, the simplest are 

5d = , 5ipi = 5a . (35) 

The first of these states that the DO-brane voltage must remain constant across the horizon 
after the perturbation. The second states that temperature must remain constant. 

Clearly, for any point r > r , Cauchy data for the perturbation equations consists of a 
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value for the ten quantities in (|34|). But because of the reduction of order at r = r , one 
cannot simply specify Cauchy data right at the horizon! Instead, one takes the five relations 
on the quantities (|34|) . obtained by setting r = r , and imposes them at r = r + e for some 
small e: these conditions are the stand-ins for the true horizon boundary conditions, and I 
believe they become equivalent in the limit e — > 0. 

To specify Cauchy data at r = r + e, one must set values for five of the quantities 
(J34|) in such a way that the other five can be determined from the boundary conditions 
described in the previous paragraph. One appropriate choice of these five quantities is 
5a, 5$, 5ip' x , 5(f2, S<f3. Of course, one must also specify r , e, and k. As usual with classical 
equations, there is a single scaling symmetry corresponding to rescaling the Planck length. 
This allows one to set tq = 1. Because the perturbation equations are linear, there is 
another scaling symmetry: qi(r) — > Ag,(r) for arbitrary A. Let us use this symmetry to set 
5a (tq + e) = 1. Let us also regard e as a fixed quantity. Then the quantities to be varied are 



where the functions are all evaluated at r = r + e. 

The strategy to find a normalizable perturbation is to vary the five quantities (jHUj) . inte- 
grate the perturbation equations, and check normalizability at some large r. The equations at 
large r have five normalizable and five non-normalizable solutions, so with the five variables 
(I36j) in hand one has just enough freedom to get rid of the non-normalizable solutions. 

The computations summarized in this section are entirely analytical, and they are in 
large part available in the form of a Mathematica notebook ^5] • 

5.3 Results of numerics 

To recap the previous section: having set r = 1 and fixed a small e, one wants to vary the five 
quantities (|36|) so as to get a solution to the perturbation equations which is normalizable for 
large r. In practice, once Cauchy data is specified, I integrated the perturbation equations 
out to a certain radius ry, and then the test of normalizability was to evaluate the sum of 
the L 2 norms squared of the six functions qi(r) over the rightmost fifth of the integration 
region of the perturbation equations, minus the rightmost hundredth — that is, roughly for 
r G (0.8rj, 0.99rj) when rf 3> 1. Let us call this region /. Then the quantity to be minimized 
is 



Such an integrated test avoids an unpleasant possibility that arises if one simply makes 
Qi{ r f) = the requirement at large r: zeroes of solutions to the perturbation equations can 
arise by making several growing exponentials cancel. Also, it's awkward to implement a test 
condition of the form qi{rf) = because there are six functions and only five variables (J36j) 
to adjust. 



k, <5$, 5(p', <5y? 2 , 5ip 3 , 



(36) 




(37) 
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One is left with the problem of minimizing the everywhere-positive function T, which 
must be numerically computed for any choice of the five variables (|36j) . If the CSC is correct, 
the absolute minimum of T is a lot smaller when (a, 9) are in the thermodynamically unstable 
region. But it's hard to find the absolute minimum of a function of five real variables which 
is computationally expensive to evaluate. I employed the Mathematica function Minimize, 
which combines several standard approaches to minimization. It samples on the order of 10 4 
points, which is actually rather few when one considers that 6 5 = 7776. I used e = 1/20 and 
rf = 5 and considered 16 different points in the (a, 9) plane, shown in figure 

In performing the numerical integrations, it was noted that when k = 0, there are static 
perturbations of the black hole which are regular at infinity but not normalizable: the 
qi(r) approach constant values. These perturbations are nearby supergravity solutions with 
slightly different charges or mass. This creates a difficulty in determining the correct normal- 
izable solutions: the function T described in ()37|) has a relatively broad, shallow minimum 
around k = 0, and numerical minimization of T is usually drawn to it even when the absolute 
minimum is at k ^ 0. One solution was to tell Minimize never to try a value of k less than 
some arbitrary lower limit (for example, 0.3). Another was to minimize not T itself, but 
rather T = e~ 1AkTf T. The choice of prefactor was arranged specially to eliminate the broad 
shallow minimum without causing T to vanish exponentially at large k: indeed, T still has 
weak exponential growth at large k. This second method proved quite effective at finding 
the normalizable solution when the CSC says there should be one. However, for reasons I 
do not understand, it seemed less good than unrestricted minimization of T itself at finding 
the least divergent solution when there is no normalizable mode. Perhaps we should not be 
too surprised: if numerics are optimized to do one thing well, they may do another thing 
less well. 

We can also understand at this point why it is difficult (though not impossible) to test 
the hypothesis that the wavelength of the stationary perturbation diverges as one passes to 
the NCFT limit. Diverging wavelength means small k, and for small k it becomes harder 
to distinguish between normalizable solutions and solutions which are regular at infinity but 
not normalizable. 

An overview of the results of numerics is presented in figure 121 Figure E3 shows a plot 
of the six scalars involved in the perturbation for a particular choice of a and 9 where an 
instability exists. The numerical results described represent on the order of 100 hours of 
CPU time on a 2.4 GHz Xeon processor (that is, a fast PC by 2004 standards). 

6 Conclusions 

The D2-D0 bound state obviously must exhibit a Gregory-Laflamme instability sufficiently 
far from extremality. The correlated stability conjecture predicts a particular threshhold 
at which the GL instability sets in: it is where local thermodynamic stability, defined in 
terms of the Hessian matrix of susceptibilities, is lost. This threshhold is described by a 
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Figure 2: Left: the minimum of T for various values of a and 9. In the final column, "False" 
means that a normalizable solution should exist according to the CSC, and "True" means 
that it shouldn't. Right: the position of these 16 points in the Q2/M, Qo/M plane. A red 
dot indicates that a normalizable mode was found, and a green dot indicates that it wasn't. 
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Figure 3: The normalizable mode for the point labeled 12 in figure El 
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curve in the plane of Q2/M and Qq/M, shown in figured Q2/M is the mass fraction of the 
D2-branes; Qq/M is the mass fraction of the DO-branes; and the remaining mass fraction is 
from the non-extremality. 

The stability threshold predicted by the CSC passes from a finitely non-extremal D2- 
brane with no DO-branes to a near-extremal bound state whose mass comes mostly from DO- 
branes. The dynamics of this near-extremal bound state is supposed to be non-commutative 
field theory at finite temperature, but the Gregory-Laflamme instability occurs for a tem- 
perature that is low compared to the intrinsic scale l/\fd of the NCFT. I conjecture that 
the wavelength of this instability becomes very long in the NCFT limit, avoiding a con- 
tradiction with decoupling arguments. But it is worth noting that for a specified D2-D0 
configuration that is close to the NCFT limit, the supergravity approximation (according 
to the CSC) predicts the existence of an instability at very low temperatures that makes 
the non-commutativity parameter depend on space (and probably time). This instability is 
difficult to see in a weakly coupled description. 

To verify the predictions of the CSC, we have looked numerically for a normalizable, 
stationary perturbation of the D2-D0 bound state for various values of Q2/M and Qq/M. It 
helped first to reduce the problem to two dimensions parametrized by r and x 1 (the direction 
in which the inhomogeneity develops). Two points in the numerical analysis deserve some 
scrutiny: 

1. Boundary conditions at the horizon were implemented in a slightly ad hoc fashion, 
motivated by the realization that the horizon is a regular singular point of the equations 
of motion when expressed in terms of a standard Schwarzschild-like radial variable. But 
the simplest of these boundary conditions are physically correct, and as a whole they 
clearly pick out finite perturbations at the horizon; so I have little doubt that the 
boundary conditions are indeed correct. 

2. The numerics were "trained" to avoid a broad shallow local minimum for small wave- 
number so as to more efficiently find the global minimum (when an instability exists) 
at finite wave-number. It is a bit disappointing that optimizing the numerical routines 
in this way made them worse at finding the least divergent perturbation in the regime 
where the CSC predicts no instability. However, it would be all but inconceivable to 
find that an instability exists when the CSC says it shouldn't: such a situation would 
amount to a violation of the Second Law of Thermodynamics as formulated in classical 
supergravity. So this optimization issue doesn't loom very large in the final analysis. 

The bottom line from the numerical study is that the predictions of the CSC are spectacularly 
confirmed, with a clean separation of stable and unstable points and scalar profiles which 
exhibit the expected exponential tails as well as some interesting structure near the horizon 
(figure EI). 

A variety of other brane bound states can be studied by similar methods. It is fairly 
straightforward to make predictions of the stability boundary using the CSC, and a typical 
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behavior is for the stability boundary to pass into the near-extremal regime in a limit where 
most of the mass comes from branes that are delocalized in some spatial dimension. Only, for 
D4-D0, where the binding energy is zero, one naturally expects a GL instability for all values 
of DO-brane charge and non-extremality. One may also consider adding angular momentum: 
it is an additional thermodynamic quantity which can be locally redistributed, so it can 
participate in GL instabilities. 2 Studying the linearized equations of motions around brane 
bound states is in general somewhat complicated, but we are fairly confident that the CSC 
would be confirmed by such analysis in all the cases we have described here. 3 I hope to 
report on these issues in the future; work on some of them is already under way |19j . 
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